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Abstract 

The aim of this paper is to apply direct methods to the study of 
integrals that appear naturally in Statistical Mechanics and Euclidean 
Field Theory. We provide weighted estimates leading to the exponential 
decay of the two-point correlation functions for certain classical convex 
unbounded models. The methods involve the study of the solutions of 
the Witten Laplacian equations associated with the Hamiltonian of the 
system. 



1 Introduction 

In these notes, we study partial differential equation techniques for problems 
coming from equilibrium Statistical Mechanics and Euclidean Field theory. In 
the context of classical equilibrium Statistical Mechanics, one is interested in a 
natural mathematical description of an equilibrium state of a physical system 
which consists of a very large number of interacting components. Consider, for 
example a piece of ferromagnetic metal (like iron, cobalt, or nickel) in thermal 
equilibrium. The piece consists of a very large number of atoms which are lo- 
cated at the sites of a crystal lattice A. Each atom shows a magnetic moment 
which can be visualized as a vector in M^. This magnetic moment is called the 
spin of the atom and represents the orientation of the atom in the lattice. The 
set S of all possible orientations of the spins, is called the state space of the 
system. Each element j of A is called a (lattice) site. A particular configuration 
of the total system will be described by an element x = ixi)^^j^ of the product 
space il, — S^. This set fi is called the configuration space. 
The physical system considered above is characterized by a sharp contrast: the 
microscopic structure is enormously complex, and any measurement of micro- 
scopic quantities is subject to Statistical fluctuations. The macroscopic behav- 
ior, however, can be described by means of a few parameters such as magne- 
tization and temperature, and macroscopic measurement leads to apparently 
deterministic results. This contrast between the microscopic and the macro- 
scopic level is the starting point of Classical Statistical Mechanics as developed 
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by Maxwell, Boltzman, and Gibbs. Their basic idea may be summarized as fol- 
lows: The microscopic complexity may be overcome by a statistical approach, 
and the macroscopic determinism then may be regarded as a consequence of a 
suitable law of large numbers. According to this philosophy, it is not adequate 
to describe the state of the system by a particular element x of the configuration 
space n. The system's state should rather be described by a family of 5— valued 
random variables or (if we pass to the joint distribution of these random vari- 
ables), by a probability measure ji on Q. consistent with the available partial 
knowledge of the system. In particular, ^, should take account of the a priory 
assumption that the system is in thermal equilibrium. 

Which kind of probability measure on is suitable to describe a physical sys- 
tem in equilibrium? The term equilibrium clearly refers to the notion of forces 
and energies that act on the system. Thus one needs to define a Hamiltonian $ 
which assigns to each configuration x a potential energy In the physical 

system above, the essential contribution to the potential energy comes from the 
interaction of the microscopic components of the system and a possible exter- 
nal force. As soon as a Hamiltonian $ have been specified, the answer to the 
question is generally believed to be the probability measure 



Here d\ refers to a suitable a priory measure (for example the counting measure 
if O is finite), /3 is a positive number which is proportional to the inverse of 
the absolute temperature and Z > is a normalization constant. The above 
measure fi is called the Boltzmann- Gibbs distribution. 

As we have mentioned above the number of atoms in a ferromagnet is extremely 
large. Consequently, the set A in our mathematical model should be very large. 
According to a standard rule of a mathematical thinking, the intrinsic prop- 
erties of large objects can be made manifest by performing suitable limiting 
procedures. It is therefore a common practice in Statistical Physics to pass to 
the infinite volume limit |A| oo. (This limit is also referred to as the ther- 
modynamic limit). The Boltzmann-Gibbs distribution does not admit a direct 
extension to infinite systems. However, when dealing with infinite systems, we 
can still look at finite subsystems provided the rest is held fixed. Indeed, start- 
ing with an interacting potential <p we can define for each finite subsystem A a 
Hamiltonian <f>^ which includes the interactions of A with its fixed environment. 
The methods for investigating phase transition for certain physical systems took 
an interesting direction when powerful and sophisticated PDE techniques are 
introduced in the mathematical technology. The methods are generally based 
on the analysis of suitable differential operators 



d^{x) 




and 



Wi') = -A+ 



— - — h Hess$. 
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These are in some sense, deformations of the standard Laplace Behrami opera- 
tor. These operators, commonly called Witten Laplacians, were first introduced 
by Edward Witten [18] in 1982 in the context of Morse theory for the study 
of topological invariants of compact Riemannian manifolds. In 1994. Bernard 
Helffer and Johannes Sjostrand [8] introduced two elliptic differential operators. 

:= - A + V$ • V 

and 

a'^^^ := - a + V$ • V + Hess$. 

These later operators, provide direct methods for the study of integrals and 
operators in high dimensions of the type that appear in statistical mechanics 
and euclidean field theory. In 1996, J. Sjostrand [13] observed that these so 
called Helffer-Sjostrand operators are in fact equivalent to Witten's Laplacians. 
Since then, there have been significant advances in the use of these Laplacians to 
study the thermodynamic behavior of quantities related to the Gibbs measure 
Z~^e~'^dx. As a simple illustration, if one is interested in the study of the mean 
value {g)j^ where 



and 




for a suitable smooth function g, one can first solve the equation 

Vg = (- A + V$ • V) V + Hess$v, 
for a C°°— solution v where the operator 

-A + V$- V 



acts diagonally on each component of v. Under suitable assumptions on the 
Hamiltonian one can see that v is also a solution of the system 

g =< g >A +v • V$ — divv. 

If it turns out that g{0) = and is a critical point of then 

< g >A= divv{0). 

Thus, the study of the thermodynamic properties of the mean value is then 
reduced to estimating the derivatives of the solution v. 

One of the most striking results is an exact formula for the covariance of two 
functions in terms of the Witten Laplacian on one forms, leading to sophisticated 
methods for estimating the correlation functions. This formula is in some sense 
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a stronger and more flexible version of the Brascamp-Lieb inequality [1]. The 
formula may be written as follow: 

cov(ff, /i) = y (aW" Vff • V/i) e-^^-Ux. (1) 

To understand the idea behind this formula, let us denote by (/) the mean value 
of / with respect to the measure 

the covariance of two functions / and g is defined by 

cov{g,h) = {{g-{g)){h^m. (2) 

If one wants to have an expression of the covariance in the form 

cov{g, h) = (V/i • w)^2(R„_Rn.e-*d^) ' (3) 

for a suitable vector field w, we get, after observing that V/i = V(/i — {h)), and 
integrating by parts: 



cov 



r{g^ h) = j {h - (/i))(V$ - V) • we-*(^)da;. (4) 
This leads to the question of solving the equation 

9-{9) = (V$ - V) • w. (5) 
Now trying to solve this above equation with w = Vu, we obtain the equation 



g-{g)=Afu 
{u) = 0. 



(6) 



Assuming for now the existence of a smooth solution, we get by differentiation 
of this above equation 

Vg = A^^^Vu (7) 

and the formula is now easy to see. 

New methods that are purely based on spectral analysis have been recently 
developed by Helffer-Bodineau [2], Sjostrand-Bach-Jecko [26]. In these papers, 
the authors studied a certain class of unbounded spin models by means of the 
spectra of the Witten Laplacian. In [26], the asymptotics of the two point 
correlation function to leading order in (3~^ was obtained under under weaker 
assumptions on the Hamiltonian. In 2003, V. Bach and J. S. MoUer [27] pro- 
posed a refined version of the results in [26] by introducing a new twisted Witten 
Laplacian to relax the convexity assumptions. 

We attempt in this paper, to study weighted estimates that lead to the 
exponential decay of the two-point correlation functions for certain convex 
unbounded systems. We removed limitations of earlier work of Helffer and 
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Sjostrand [8]. They only treated the one dimensional case {d = 1) under the 
artificial restrictions 

||Hess$(a;)||£(,<.) < C 

and 

||Hess$(a;) - I||£(ioo) <S<1, 
for all weight function p on Z/mZ satisfying 

g-« ^ ^i!±ll < for some k>0. 

These conditions are too restrictive for many important applications, while my 

conditions arc considerably more flexible. In particular, the conditions in my 
work are suitable for treating the d-dimensional nearest neighbor Kac model, 
where the potential is given by 

^i^) = — 2 ^ In cosh 

for v > smaller than some value Uo to be determined. 

In section 2, we give a motivational background on the origin of Witten's 
Laplacians. 

In section 3, we give an outline of the operators and equations involved in 
the Witten Laplacian method. 

In section 4, we discuss preliminary results on Hilbert space methods for 
elliptic PDE's. 

In section 5, we provide a rigorous discussion based on Hilbert space methods 
for the solvability of the corresponding Witten Laplacian equations. 

In section 6, we illustrate the family of Hamiltonians discussed in section 3 
and 5 through an example of the type introduced by Marc Kac [20]. Section 
7 is devoted to the study of the exponential decay of the two-point correlation 
functions for models of Kac type in the convex case. We shall establish weighted 
estimates leading to the exponential decay of the two-point correlation functions. 

In section 8, we shall apply our methods to the d-dimensional nearest neigh- 
bor Kac model. 
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2 The Witten's Laplacians 

In 1982, Edward Witten published an article [18] on Supersymmetry and Morse 
theory relating invariants of a Riemannian manifold M with some indices of a 
Morse function <I> € C°°(M). For this, he introduced the Witten derivative d$ 
and the Witten coderivative dj by simply setting 

d$=e~^de^ and dj=e?d*e~^ , (8) 
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where d and d* are the exterior derivative and exterior codcrivative respectively. 
The Witten Laplacian is then defined to be the associated second order operator 

W$ = (d$+d|)' (9) 
= d$di + djd<i, (10) 

acting on the exterior algebra bundle of the cotangent bundle of M as the 
standard Laplacian does. 

Choosing a local orthonormal frame field ei, and denoting by e^, e*^ 
its dual cofranie field, d and d* could be easily represented in terms of the 
Riemannian connection V as 

d = e*AVe, and d* = -'i^.Ve,. (11) 

where denote the interior product with respect to Cj (sec [58] for more 
details). Here and in the rest of this section, we use the Einstein summation 
convention namely, an index occurring twice in a product is to be summed from 
1 up to the space dimension. We consequently have 

d$ = e'AVe.+e^^ and = -i(e^.) V,^. + i(e^-)^ (12) 

where ^-i^i^... denote the components of multiple covariant differentiation rela- 
tive to the local frame field ei, ...,ef;. 

= Ve,Ve.$-Vv..e..$ (13) 

Since e' AVp and i(e -) V„ do not depend on the choice of the local orthonormal 
frame and coframe field we may assume that ei,...,ed comes from a normal 
coordinate centered at an arbitrary point, and consequently have 

Ve^e'A = Ve,i(e^.) = 0. (14) 

Now using (10), (11), (14) and the fact that 

e' Ai{ej) + i{ej)e'A = 6ij, (15) 

we have 

A+^i^+^(e^Ai(e,) - i(e,)eV). (16) 

In the case of R" where covariant differentiation becomes standard differentia- 
tion, the Witten Laplacian on 0-forms acting on a smooth function / gives 

Wi°)/ = -Af+^^f-^f (17) 
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The Witten Laplacian on one-forms acting on a one form 

u = vf'{x)dx'^ 

gives 

W^^\= Au+^^^u-^^u + 2^^dx'M^u. (19) 
4 2 2 a^.^ 

Identiiying one-forms with vector fields in M" (1-12) becomes 

|2 



W^^^u = j^- A+^^L - O u + Hess$u. (20) 

The tensor notation simply means that the operator —AH ^ — acts 

diagonally on each component of the vector field u. Let us also point out that 
the identification between forms and vector fields is a common practice in Rie- 
maniann geometry and is done via the metric tensor. 

As first observed in [8] by Bernard Helffer and Johannes Sjostrand, these 
Laplacians provide new methods for solving problems coming from Statistical 
Mechanics. The methods are generally based on the analysis of the differential 
operators 

:= - A + V$ • V (21) 

and 

A^i) := vl^ ^ 7rf+Hess$. (22) 

These two elliptic differential operators for which a Predholm theory can be 
developed are equivalent, as observed in [13], to Witten's Laplacians W^^ and 
wj,^^ respectively where 

V$|^ A$ 



and 



W^^^= -A+^ (23) 



Wi'' = ( - A+^^^ - ^ I ® I + Hess$. (24) 



Indeed, it only suffices to observe that 

iy^)=e-*/2o4)oe*/2 (25) 

and the map 

u I — > e 2 u 

is unitary 
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3 The Basic Equation 



For any finite domain A of Z'^, we shall consider a Hamiltonian (f>A of the 
phase space M^, satisfying conditions that will guaranty the solvability of the 
corresponding Witten Laplacian equations. We shall also consider a slowly 
growing source term g, to ensure that the solutions have suitable asymptotic 
behavior. 

We shall first establish the solvability of the equation 

^^'^L^M (26) 

by means of Hilbert space methods. The method consists of determining an 
appropriate function space and an operator which is a natural realization of the 
problem. In this particular problem, the function spaces to be considered are 
the Sobolev spaces B^{R-^) defined by 

B|(M^) = {ue L^{R^) : Zid"u e L'^{R^) V ^ + |a| < fc} . 

where 




|V*| 



(27) 



2 

These are subspaces of the well known Sobolev spaces W'''^{M.^), A: e N. 

The vital tool in the Hilbert space approach to elliptic boundary value prob- 
lems is the celebrated Lax-Milgram theorem. The essence of the method is the 
interpretation of the problem in a variational sense involving a bilinear form 
defined in a natural way by the problem and acting on the appropriately chosen 
function spaces. 

In general, the Hilbert space method for elliptic differential equations uses 
the Compact embedding theorem for Sobolev spaces. This is a fundamental 
step in the method in order to be able to apply the Predholm alternative. Since, 
in the context of our problem we are dealing with unbounded domains, the 
classical results regarding the compactness of the embedding 

W'''P{n, dx) ^ LP{n, dx) (28) 

for suitable O are no longer valid. However, In the case where the spaces are 
taken with respect to the weighted measure e~*dx, with a suitable we have 
the following result due to J-M. Kneib and F. Mignot [ll]lem.5. 

Lemma 1 // $ satisfies the condition 

36* e (0,1): lim 6l|V$(a;)|^-A$ = oo 

|a:|— >oo 



then 

is compact. 
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Here and in the sequel, rf/x will denote the Gibbs measure 

djj, = Z~^e~^dx, 

Z = [ e-^dx. 

and H^(ij) denotes the weighted Sobolev space 

H''{^x) = {u€ L^{M.^, dn) : d"u G L'^{R^, d^i) V |a| < k) . 

Proof. We shall prove that every bounded sequence in has a convergent 

subsequence in Let {uk} C = ^^(1^"^,^^) be such that 



llufell^i < vM for every k and some M > 0. 

For any i? > 0, denote by -6(0, R) the open ball centered at with radius R. It 
is clear that H^{R^, dfi) C H^{B{0, R), djj). Hence {uk} is a bounded sequence 
^"1(5(0, R), rf/x). Moreover 



/ u1.dx + I \Duk\dx 
Jb(o,r) Jb(o,r) 

Jb(o,} 



lB{0,R) JB{0,R) 
< C<E, 



\Duk\e "^dx 

B{0,R) JB(0,R) 

This implies that is a bounded sequence in H^(B{0, R)). Now using the 

standard Sobolev compactness embedding theorem for bounded domains with 
nice boundary (see [3]), we get the compactness of the embedding 

H^{B{i},R)) ^ L^{B{0,R)). 

Therefore, one can find a subsequence {wfe^ } of {ufc} such that Ui-^ converges in 
L2(B(0,i?)).Wc shall prove that {wfej is Cauchy in L'^{R^,dn). Let r] > O.The 
assumption of the lemma implies that 

(:=|V$|^-(1 + 77)A$ (29) 

is positive in a neighborhood of oo when 6 = (l + ?7)~^. 

[ \ukj-Ukif e~'^dx < [ \ukj-Ukif e^'^ + [ /' ^''^ e'^SS) 



'\x\<R «'|x|>iJ , „,„ „, 

' ' ' \ B{0,R) 

2 /' C\ukj-Uk, 



J\x\<R J\t 



|2 



To estimate the last term of the right hand side of this last above inequality, let 
s > and choose R large enough so that 

■^f ^ > 4M(2 + >7 + r?-i) _ 

M'^ \ B{0,R) ~ e 
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Now introduce the vector fields 

Xj = dj (32) 

and their formal adjoint in L^{^) 

= -9, (33) 

one has when u G C^(K^) for their sum and commutator 

{Xj + X;) u = 'l>,^u (34) 

and 

[X„x;]u = ^,^,.u. (35) 
It is then straightforward to see that 

{[X„X;]u,u)^ = \\x;u\\l-\\X,u\\l (36) 

\\{Xj+X;)u\\l<(^l + ^y\X,u\\l + {l + n)\\X*u\\l , ye>0 (37) 
so that a linear combination of these formulae gives for any ry > 

((I V$|2 - (1 + 7j)A<^)u,u)^ < (2 + 77+7,-1) (\\X,u\\l + - + \\Xmu\\l) . (38) 
Thus, 

{iCu,u)^<{2 + rj + rj-')\\u\\l,^^) (39) 

Because Cj^(M'^) is dense in H^{ii), this inequality is valid for all u € H^{fi). 
Now applying (39) with u replaced by Ukj-Uk,, (31) gives 

V ' ' J|x|<fi' ' 4M(2 + j? + r?-i) 

< C$ / |ufe^._UfeJ^ +£ 

^|x|<iJ 

The result follows from the convergence of the subsequence {ukj } in L^{B{0, R)). 
■ 

The Lemma above; indicates the direction towards the assumptions needed 
for the Hamiltonian $ = $a- 
Assumptions on $. 

Recall that A is a finite domain in Z'' . We shall assume that $(x) e C°°(M^) 
satisfying: 

1. lim I V$(x)| = oo 

2. For some M, any 9"$ with |a| = M is bounded on K^. 

3. For \a\ > 1, there are constants such that \d"^{x)\ <Ca(l + | V$(a;)|^ 

4. Hess$ > 6 for some < 5 < 1 
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4 Preliminary Results on Hilbert Space Meth- 
ods For Elliptic PDE 

A bilinear form with domain a complex Hilbert space, is a complex-valued 
function a defined on Hx H which is such that a{u, v) is linear in u and conjugate 
linear in u.Tlie inner product (•, on H is clearly a bilinear form; we shall 
denote it by 1 (•, •) .The form a + XI will simply be denoted by a + A 

(a + A) (u, v) = a{u, v) + X{u, v)h- 

The adjoint a* of a is defined by 

a* (it, v) = a{v, u) 

and a is said to be symmetric if a = a* , i.e. for all u,v G H 

a*{u, v) = a{v, u) = a{u, v). 

A bilinear form is said to be bounded on if x iJ if there exists a constant M > 
such that 

\a{u,v)\ < M \\u\\^j \\v\\jj for all u,v&H. 

A bilinear form a is said to be coercive on H if there exists a positive constant 
m > such that 

|a(u, m)| > TO ||w,||^ for all M, w G 

We shall say that a Banach space W is continuously embedded in a Banach 
space X if there is a bounded operator E : W ^ X which is one-to-one. We 

call E an embedding operator. We shall say that W is densely embedded in X 
if R{E), the range of E is dense in X; and we shall write 

W -^l X. 

If X is a Banach Space, the set of all linear conjugate functionals on X shall be 
denoted by X* and is called the conjugate space of X* . 
Suppose X, Y, W, Z are Banach spaces such that 

W^gX and Y^gZ*. 

Let a{w,z) be a bounded bilinear form on W x Z. We can define two linear 
operators connected with a{w,z). The first which we shall denote by A, is an 
operator from X to Y. We say that x G D{A), the domain of A and Ax = y if 
X e R{E), y €Y and 

a{E~'^x, z) = Fy{z), for all z € Z. 

Since R{F) is dense in Z* , the operator A is well defined. We call A the operator 
associated with the bilinear form a{u,v). 
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The second operator, which wc denote by A, is from W to Z* . We define it as 
follows. Fix w G W, a{w, •) G Z* , it is bounded because the bilinear form a is 
bounded. We define Aw to be a{w, •). A is clearly well defined and will be called 
the extended linear operator associated with the bihnear form a{u, v). It can be 
shown that A and A are related in the following way: 

A = F-^AE-^ 

The fundamental tool to investigate the operator A is the Lax-Milgram theorem 

Theorem 2 (Lax Milgram) Let a be a bounded coercive form on a Hilbert 
space with bounds m and M as above. Then for any F e H* , the adjoint of 
Ho, there exists an u € Ho such that 



The map A:u\ 



a{u, v) = {F, v) for all v G Hq 

F defined above is a linear bijection of Ho onto H* and 



m < 



A 



< M, 



M- 



< 



< m 



Proof, (see [3]). 



Corollary 3 For any choice of F € H* there is a unique vector u G Ho satis- 
fying 

moreover, the isomorphism A~^ from H* onto Hg defined by A~^F = u verifies 



A-^F 



\F\ 



Ha 



Next, we apply the Lax-Milgram theorem to the situation where the Hilbert 
space Ho is continuously and densely embedded in another Hilbert H. 

Lemma 4 If H is a Hilbert space and W is a Banach space continuously and 
densely embedded in H with embedding operator E, then H can be continuously 
and densely embedded in W* with embedding operator F satisfying 



{x, Ew)jj = Fx{w) 



X € H, and w €W. 

{x,Ew)fj is a conjugate linear 



Proof. For each x G H, the function x* : w < 
functional on W and 

\x*{x)\ < \\x\\h\\E\\\\w\\^. 

Hence x* G W*. Define the operator F from H to W* by Fx — x*. Clearly, F 
is linear and bounded. It is also one-to-one since R{E) is dense in H. Finally, 
suppose x*{w) = for all x* G RiF). Then {x,Ew)^ = for all x e H. Thus 
Ew = and consequently w = O.This shows that R{F) is dense ■ 

Now let the Hilbert space Ho be continuously and densely embedded into 
another Hilbert space H with embedding operator E. By the lemma above, 
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H can be continuously and densely embedded in H* with embedding operator 
F.We obtain the scheme 

which is referred to by saying that {Hg, H, H*) is a Hilbert triplet. Notice also 
that if the embedding E is compact, then so is the embedding 

Ho --^^ K. 

Returning to the biUnear form on Ho, we weaken the notion of coerciveness as 
follows: Wc say that a bilinear form a{u,v) on Hg is coercive relative to H, if 
there exists some A > such that ax{u, v) = a{u, v) + X {u, v)fj is coercive, i.e. 

a{u, u) + X \\u\ffj > Uo \\u\\% for u G Ho and some Uo > 0. 

If this last inequality above holds, then by Lax-Milgram, the extended linear 

operator Ax associated with the bilinear form ax{u,v) has a bounded inverse 
: H* Ho, moreover Axu = Au + XBu, where A is the extended operator 
associated with the bilinear form a{u, v) and B the extended operator associated 
with the inner product {u,v)^ . 

Now Let q G H* and consider the equation 

u e Ho, Au = q (40) 

(1.8) can now be written as 

uG Ho, u- XA^'^Bu = z (41) 

with z = A'^^q.We now claim that the compactness of the embedding E implies 
that of the operator A^^B : Ho ^ Ho is compact. Indeed this follows from 
the fact that B is bounded and A^^ : H* Ho is compact. By the Fredholm 
alternative (see theorem 3 below), (1.9) is uniquely solvable for any choice of 
z G Ho'ii and only if m = is the unique vector of Ho satisfying u—XA^^Bu = 0. 
When this is the case, the linear operator z i — > u defined by (1.9) is bounded 
from Ho to Ho- Summing up, we have the following theorem 

Theorem 5 Let {Ho, H, H*) be a, Hilbert triplet with Ho compactly embedded 
in H, let a{u,v) be a bounded bilinear form on Ho coercive relative to H. Then 

u G Ho, a{u, v) = q{v) for v G Ho 

admits a unique solution u for any choice of q G H* if and only if it admits a 
unique solution u = for q = in which case the solution u satisfies 

Mh^<CMh^ 

with C dependent only on A. 
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Theorem 6 (Predholm Alternative) Let T be a compact linear operator on 
a Hubert space V and consider the equations 

ueV, u-Tu = f (42) 

veV, V* - T*v* = g (43) 

where T* the adjoint operator of T. Then the following alternative holds: 

(i) either there exists a unique solution of (42) and (43) for any f and g in V, 

or 

(a) the homogeneous equation 

u-Tu = 

has nontrivial solutions. In that case the dimension of the null space of I — T 
is finite and equals the dimension of the null space Af* of I — T* Furthermore 
(42) and (43) have solutions (not unique ) if and only if 

{f,v*) = 0, \/v&M* 

and 

{g, v) =0, VveAf 
Af being the null space of I — T 

Proof, .see Yosida [17], (X - §5) ■ 

Remark 7 Assumption 2 implies that ^ is a slowly increasing function. This 
assumption is made to rule out any possibility of exponential growth for 

5 Solvability and Regularity of The Basic Equa- 
tion 

Theorem 8 Let K be a finite domain in . If ^ satisfies assumptions 1-4 
above, then for any C°°— function g satisfying 

\D"g\<Ca{l + Z^)<'- (44) 

where 

7 |V$| 



2 ' 

a e NI^I with some Ca and some qa > 0, there exists a unique C°°— function u 
solution of 

{ Afv = g-{g) , , 

{ * ^ ^^'l^m (45) 
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Proof. (Existence) Wc shall work in the unweighted space L^(K^) and with 
the Witten-Laplacians ensuing after the unitary transformation. 
Under the unitary transformation, 

A''^\ = g-/g\ inlRA 

is equivalent to 

W^°^u = q in 

where 

u = e-^/\ and q = e-^/'ig - {g) ) G L\r''). 

Let 

B|(rA) = {u€ L'^iR^) : Z^^d'^u G \f I + \a\ < k} 

(here d"u is taken in the distributional sense in M.^). 

Denote by ^.(M^) be the closure of Co"(M^) in ^^(M^), and let b be the 
bilinear form on ^.(M^) defined by 

b : Bl^^iR^) X Bl^^iR^) ^ M 

with 



i{u,w) = / 
Jr 



Du ■ Dwdx + / I - — ^ — J uwdx. 



Because wc have in mind to apply theorems 6 and 7 above, we need to check 
boundedness and coerciveness of b. 
Boundedness: After observing that 

A$ < C(l + |V$|^)i/2 

< C(l + I Vc&l^), 

it then follows immediately from Cauchy-Schwartz inequality that 

\h{u,w)\ < Q;o||w|Ib1(irA) IkllBi(lRA) 

for some constant Uo > 0. 
Coerciveness: 

[ \Dufdx = h{u,u)- I I ^-^P- - ^] \uf dx 

f 2 f 2 f 2 

/ \Du\ dx + |Z$w| dx = h{u,u)+ / — |m| dx 
Jra Jra Jra 2 

f (A^f 2 1 /■ 2 

< h{u,u) + e / — - — \u\ dx + — \u\ dx 



< b(M,u) + C£ / \Z^u\^ dx + ( Ce + —] I \u\^dx 
RA V 4ey Jra 
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choosing e such that Ce < 1 and adding J^a |u| dx on both side of this above 
inequahty, we immediately get 

^ ll"lls4(MA) < b(u, u) + 7 ||u||i2(KA) (46) 

for some positive constants 5 and 7. 

This shows that the bilinear form h{u,v) is bounded and coercive relative 
to L2(kA)_ 

Observe that B\ ^(M^) is densely embedded into L^(M^). Now considering 
the Hilbert triplet 

(i?i_^(M^),L2(M^),i3-;(M^)), (47) 

where -B~^(M^) denote the conjugate space of ^(M^). 
We need to check that the embedding 

is compact. This follows from Lemma 1 by simply observing that 

<*(M^) C U^' {H\n)) 

and the fact that f/$ is a unitary operator. 

Let be the bilinear form in ^(M^) defined by 

B^{U, W) = h{U, W) + 7 {U, W) j;,2(RA) 

and 

be the extended linear operator associated with the bilinear form B^{u,w).We 
have 

A^u = Au + ')Bu, (48) 

where A and B are the bounded bilinear forms associated with b and (•,-)i,2 

respectively. 

Note that the equation 

u e B^^iM.^) Au = q 

is the variational interpretation of the equation 

W^°^u = q in K^. 

By theorem 1 (Lax-Milgram) , the boundedness of B..,, and the coercivity condi- 
tion 

B^(u, u) > 5 ||m||bi(kA) Vw G Bl^{R^) 
guarantee that A^ has a bounded inverse 

: B-liR^) - Bl^^iR^). 
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Now using the fact that 

A^u = Au + jBu, 

we can write the equation 

u e Bl^{R^) Au = q 

as 

u G <$(M^), u - ^A-^Bu = z (49) 

where 

z = A-\. (50) 
As in the preliminary, because the injection 

Bl^iM."^) ^ L2(KA) 

is compact, the operator 'yA~^B : i?^ ^(R^) — > S^_^(M^) is compact. Moreover, 
the boundedness of ^A~^B implies that 

(7!;^^)* = (7(5-1^)*)"' (51) 

= 7(i;(^-i)y' (52) 

= 7(i;(s*)"') (53) 
= lA-'B. (54) 

Let us also point out that the self-adjointness of A^ and B follow from the fact 
that they are both associated with symmetric bilinear forms. 
Now observe that 

ker(J - ^A-^B) c ker A. (55) 

We now claim that 

kcri= |5e-*/^<5eM}. (56) 



Indeed if Au = , then b(u, u) = 0. Hence 

2 

= 



9. . „ 

which would imply that m is a solution of the equation 

5^ + — ju = 0. 
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One can then easily see u must be a constant multiple of e~*/^.We have in mind 
to apply the second part of Theorem 7 (Fredholm alternative). This brings us 
to check orthogonality of q with ker(7 — 'yA~^B). Let i5 e M, 



1,2 (RA) 



/ <5e-*/^e-*/^(,-(5)^ ) 

Hence using part (ii) of theorem 3, we conclude that the equation 

Au = q 

is solvable therefore 



(.9), 



(57) 
(58) 

(59) 
(60) 



is solvable in the weak sense. To complete the proof of theorem 4, we need to 
prove that the L^— solution constructed above is a classical solution. ■ 

Regularity: Next, we shall prove that the weak solutions constructed above 
are actually classical solutions. The proof is based on the method of difference 
quotient. 

Theorem 9 (S'^-regularity) Given q € B^"^(M'^) for k = 0, 1,2, a solu- 
tion u G Bl,^{R^) of 

Au = q 

IS an element o/B|+^(M^) and we have the estimate 



(61) 



Au 



+ u 



B|(RA) 



(62) 



for allue Bl+^{R^). 

Proof. We first establish the result when A; = 0. We have 



-u , u 



L2(RA) 



1,2 (RA) 



^ ^ (RA) ||u|Il2(ra) 

< eC (RA) + ^ ||M||i2(RA) 



(63) 
(64) 
(65) 



Thus, for u e Bl^ 



Au, u 



\Du\ 



L2(RA) 

2 



> ||-D'U||j^2(RA) + ||-^4>m||i,2(rA) - eC (RA) 



c_ 

4e 



IrIIl2(MA) • 
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Choosing e such that eC < 1, we get 



Au,u) > C||U||^WRA^ -C\\U\\\2„^. 



Hence 



< C 

c 



< 



Au 



Au 



2 



|w|Ib1(rA) +C||W||^2 



(RA) 



+ II'"'IIb4(RA) + C ||u||i2(RA) . 



Again choosing e appropriately, (eC < 1) we finally get 



< C 



Au 



+ C\\u\\U 



Now assume that for u € ^^^^(R^), iu = g € S^-^(M^) implies u e S|+i(M'^) 
and that 

+ ll^llBimA-) . (66) 



3j-i(RA)) 



Suppose now that u e ,i,(E^), Au e S|(M^). So we know that u e S^+^(M^) 
and we want to estabUsh that u e B^'^iM.^). 
Because 



replacing u by D-'u in inequality (66) we get 



|-Dfu||j3fe+i(gA) < C 



< C 



where 



AD^iu 
DT^Au 

:= 



g,_l^g^^ + FiW||B|(RA) 



+ WuD^XJ 



Bi-i(RA) 



iv$r A$ 



Bj(RA) 



4 2 

Now letting /i ^ and using assumption 3 on $ we get 



ll-Djulljjfc+i^jgA^ < C 
it then follows that 



Au 



BJ(RA) 



IfIIb|(RA)) + IfIIsJ+1(RA) 



Awesi+i(K^). 



It then only remains to prove that Z^'^u G L^(]R^). To see this first observe 
that 

Z^u = Au + Au+ —^u. (67) 
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Here, the Laplacian is taken in the distributional sense. Multiplying by on 
both sides of this last equality, we obtain: 

Z^+\ = Z^Au + Z^Au + Z^^u. (68) 

The first term of this equality is in L'^{R^) because Au e B|(M^). That the 
second terms also belongs to L^(M^) follows from the fact that DjU e B^^^ (M^). 
Finally to see that the last term is an element of L^{M.^), we use assumption 3 
on $ to get that 

^ < (69) 

< C{^ + Z^), (70) 
and use the fact that u e B|+^(]R^). ■ 

Proposition 10 (C°°— regularity) The weak solution u ofW^^u = q is an 
element o/C°°(K^). 

The proof of this proposition use the general Sobolev inequalities theorem 
given below. 

Theorem 11 (General Sobolev Inequality) Let U be a bounded open sub- 
set o/K", with a C^ — boundary. Assume u G W'''P{U) where 

W'''P{U) := {u e Ll^iM.") : d°'u e (M") V \a\ < k} . 

If 

k> - 
P 

then u e C'=-[$]-^'^(i7), where 



1 = 



- 1 — - , if - is not an integer, 
any positive number < 1, if ^ is an integer. 



Here C'^'"{U) is the Holder space consisting of all functions u E C^{U) such 
that 



l'"llc'=.-(£7) X! sup |(9'^^u(a;)| + ^ sup 



d^u{x) — d^u{y) 



\x-y\ 



< DO 



Proof, see [3] ■ 

Proof of proposition 11. Because q G C°°(]R^) , we have u € B^^^{M.'^) Vfc, 
which implies u e H''{V) (= W^^'^{V)) Vfc and VF CC M'^. Now choose fc G N 
such that fc > |A| . Then the theorem above implies that u G C''''^{y) for some 
< 7 < 1. Consequently, u G C'^iV) for an arbitrary big enough fc and for any 
FCCM^ ■ 

Now that we have enough smoothness, we can make the following remark 
which completes the proof of theorem 9. 
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Remark 12 A simple integration by parts argument shows that u is in fact a 
strong solution. It satisfies 

W^^\ = q 

pointwise almost everywhere. Using the unitary transformation and taking gra- 
dient on both sides of 

we get 

A'^Wv = Vg. 
If q is a smooth vector field satisfying 

\d"(i\<Ca{l + Z^y'' for some qa>Q, (71) 

then one can show as above ( this time using uniqueness result of the Fredholm 
alternative) that the equation 

4^^v = q 

has a unique weak solution. 

A^^^'Vv = V(? would then imply that two solutions of 

Afv = g-{g)^^^^^ (72) 
must differ by a constant. Thus the problem 

[ Afv = g-{g)^^^^^ 

has a unique solution. This ends the proof of theorem 9. 



6 The Kac-like Model 

In this section, we propose to illustrate the results above through the study of 
a more specific family of classical unbounded spin model related to Statistical 
Mechanics and is given by 



$(a;) = $a(x) = — + ^{x), X e K^. (73) 
Here we have used the notation = x ■ x. 

The model that was originally suggested by M. Kac corresponds to when ^ is 
given by 

^{x) = —2 In cosh y ^ (^i + ^j) 

where u \s a. small positive constant. 

Other aspects of this family of potentials are studied in [8] in the one dimensional 
case. 
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Definition 13 The lattice support, Sg of a function g on is defined to be 
the smallest subset T of A for which g can be written as function of xi alone 
with I e r. For instance, if g = Xi, Sg = {i}. 

Under the assumptions 

|a"V*|<Ca, VaeN''^!, (74) 

Hess$ >S>0, 0<^<1, (75) 

One can check that <I> satisfies the assumptions 1-4 in section 3. 

Let g be a smooth function on where F is a fixed subset. We shall use 
the notation 

xi: = (a;i)ieE 

if S is a proper subset of A and shall also assume that Sg = T. Now define the 
function g on R^ by 

g{x) = g{xr), x e K^. 
If there is no ambiguity we shall identify g with g. 

Wc propose to prove that if in addition to the assumptions above on the 
functions 4* and g are compactly supported and g satisfies, 

\d°'Vg\<Ca, VaeN'^^l, 
then the solution v of the equation 

{-Av + V$ ■ Vv = g - (g) 

constructed in section 5 satisfies 

a"Vw(a;) ^ as |x| ^ oo VaGNl^^L (77) 
Recall that under a suitable change of variables, the equation 



could be written as 



4'^^^ = "^9 (78) 



^- A+-^^-^ — j (g) u + Hess$u = q (79) 

where 

u = e-^/'^Vv and q = e-*/^ (80) 

Let Bi = Br^ (0) C M^denote a large balls centered at zero with radius Ri and 
containing the support of ^ in K^. We also consider a ball B2 = Bji2{0) C RJ' 
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of radius R2 > Ri containing the support of g in M^. The support of g in MA is 
then contained in the cyhnder 

B = B2X 
In = R^\B we have 

-^4-?+^)"-" (81) 

u = cp on dB (in the trace sense). 

Here is a C°°— vector field on dB and m = |A| . 
Since the operator 

x"^ m ^ 

(82) 

acts diagonally on u, we can work component by component and the situation 
is reduced to the scalar case 

-^+T-2+T = " """^ • (83) 

u = If on dB (in the trace sense). 

Having reduced the problem to a Dirichlet type for the Schrodinger operator 

-A+^ - I + 1, (84) 

we shall need some results on the decay of cigcnfunctions of the corresponding 
Schrodinger operator. We need the following lemma: 



Lemma 14 The fundamental solution £ S iS'(M'^) of the operator —A+k^ [k > 
0) exists and is unique. It is spherically symmetric, is an element of C°° {M.^\{0}) 
and has the following asymptotics as \x\ ^ 00 : 

f(a;) = C7|a;|('^) e-'^l^l (1 + 0(1)) (85) 



In the Lemma, S'{M.^) denotes the space of tempered distributions on M.^ 
Proof. Consider the equation 

{-A + k^)£{x)=5o{x). (86) 
Taking Fourier transform, we get 

{-A'+k^)£{x) = 5^). (87) 

equivalently 

%/2 



(x2 + fc2) f (a;) = (27r)-'"/' (88) 
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which impUes 

(-0 ^-"l/2 

^(-) = W- ^''^ 

The uniqueness and spherical symmetry follow since 

£{x) = {2Tr)-"'^^£i^). (90) 

Furthermore, if .x 7^ 0, the smoothness of £{x) follows from the regularity theory 
of the elliptic equation as discussed above in section 5. 

(- A + £{x) = in R^\{0} (91) 

for a; 7^ set S{x) = f{r) where / e C°°(K+) and r = \x\. (91) becomes 

-f"{r)-—^f'{r)+k'f{r)=0. (92) 

Set f{r) = a{r)g{r). Plugging this in (92) and setting the coefBcient of g'{r) 
equal zero gives 

THj — 1 

2a' + a = 0. (93) 



Take 
Then 

and (92) takes the form 



a{r) 



r 

m— 1 



m— 1 

/(r) = r 2 g{r) 



g"{r)-k'{l + O{^))g{r)=0 (94) 

Now using classical results on the asymptotics of the solutions of the Schrodinger 
operator (see [65]), we discover that 

ff±(r) = Ce±'='-(l+o(l)). (95) 

Hence the asymptotics of the solutions of (92) are 

/±(r) = Cr-^e^'^^'il + o(l)). (96) 
Since £{x) = f{\x\) G »S'(M^), we conclude that f = f- and the result follows. 



Theorem 15 Let Q be any exterior domain in containing a neighborhood 
of infinity with smooth internal boundary. Let the potential v{x) G C°°{fl) and 
satisfy 

lim inf ■u(a;) > E (97) 

\x\—^oo 
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and let tp be a smooth solution of the problem 



p = if on dfl ^^^^ 

where X < E and ip is a smooth function on dfl. Then the following estimate 
holds: 

Hx)\ < ae-V(«-^-^)/2kl (99) 

for any e > 0. 

The proof of this theorem uses the foUowing lemma 

Lemma 16 (A Maximum principle) Let k > 0, T, an open subset o/M^, 
and u e C2(S) a function such that 

(-A + A:^) u = / < mS. (100) 

Then u cannot have a positive maximum in S. 

Proof. If e S is a maximum point and u{xo) > 0, then 

Au(a;o) < 0; (101) 

this contradicts (100). ■ 

Proof of Theorem 7. Let be a real solution of the equation 

Hip = Xtp in O. (102) 

where 

H = -A + v{x). 

We obviously have 

A{(p^) =2A(p-(p+2\Vipf (103) 
Hip = X<p gives —A(p = (A — v{x)) (p which implies 

A{ip'^) =2{X-v{x))(p'^ -2\Vipf (104) 

adding 2(5 — X)ip'^ on both sides of this equality, we obtain 

[- A+2(6 - A)] (fi^ = -2 {v{x) -b)(p^-2 \V(pf . (105) 

Choosing X < b < E the right hand side of (105) is non-positive for large 
enough. Now set 

u{x) = ip^{x) - M£{x) (106) 
where £{x) is the fundamental solution of the operator — A + fc^ with 

k = v'2(6- A). (107) 
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Choose R so large that £{x) > and v{x) > b for |a;| > R. Now choose M so 
large that u{x) < on {x £ Q : \x\ — R] . We shall prove that 

u{x) < (108) 

on {x € Q : \x\ = R} from which the theorem will follow. Substracting from 
(105) the equation 

[- A+2(6 - A)] M£{x) = 0, (109) 
we find that (100) is satisfied for ^(a;) with 

/ = -2 {v{x) - 6) _ 2 I Vv?|^ , for \x\ > R. (110) 

We then apply the maximum principle in each connected component of the 
subset 

ClH,p = {xGTi:R<\x\<p} (111) 

to the function 

u^x) = j u{x- y)r]Ay)dy (112) 

X 

where r]^{x) = e~"^r]{ — ) and t]{x) is the moUifier. Recall that rj^x) is given by 



T]{x) 

We indeed have 



e (i-l^lO if |a;| < 1 
otherwise. 



(-A + e) u'^r^ J fix - y)vAy)dy < (113) 

u e L^{M.^) implies that u'^{x) ^ as \x\ oo. Set 

Mp{s) = max \u%x)\ (114) 

since u{x) < for a; G {x £ : \x\ = R} , using the fact that u'^{x) ^ u{x) as 
£ — > on {a; e O : |a;| = i?} , we conclude that u^{x) < on {a; e O : |x| = i?} 
for small sit then follows from lemma 17 that 

u%x) < Mp{e) for x e Clu^p. (115) 

Letting p — > oo, we get 

u^{x) < for a; e O and \x\ > R. (116) 

Now since 

u^(.x) =1 u{x) as e -> (117) 
in every relatively compact subset of {x £ Q : \x\ > R} , it follows that 

u{x) < for X e {a; e O : |a;| > i?} . (118) 
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Corollary 17 If v{x) — > oo as |x| oo, then for any eigenfunction ip of the 
boundary value problem in theorem 7 satisfies, the following estimate 

Hx)\ < C^e-^l-l (119) 

where a> is arbitrary and Ca > 0. 

Theorem 18 (Helffer-Sjostrand [8]) The L'^— solution u of 

Tfl \ 

"^+T-¥ + T = ° iE) (120) 



satisfies 



where 



u = if on dB {in the trace sense). 



u{x) = e-~ \x\~^''^ h{x) (121) 



d^h{x)=0{\x\-^^^) V/JeN". (122) 

Using the change of variable v = e*/^u and applying this theorem to each com- 
ponent of u, we obtain 

Corollary 19 The L"^ — solution v of the system 

(- A + V$ • V) V + Hess$v = V(/ inM.^ (123) 

satisfies 

lim a"v(a;) =0 Va e N™. (124) 

Proof of theorem 8 (Sjostrand [14]). Denote by 

K : C°°{dB) C°° (B") (125) 

the operator that assigns each boundary value the corresponding solution. Since 
by theorem 16 

lim u{x) = 0, (126) 

the maximum principle implies that K is monotone increasing. Indeed, Kg > 
whenever g > 0. This implies that the operator K is increasing and that Kg < 
sup 5, if sup 5 > 0, Kg > g if inf g < 0. Let 

Uo = Kl (> 0) (127) 

which is a radial function i.e. 

Uo{x) = Uo{\x\); (128) 

with 



c,2 , m - 1\ m 
-dr - I ' ^'^ + T ~ 2" 



Uo(r) = 0, Uo{R) = 1. (129) 
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We perform the Liouville's transformation 



u„ = r 



i)/V(r) (130) 



to get rid of the term involving d,.. We finally get 



/(r) = 0, f{R) = 

(131) 



, !_ _ (m - 1) (to - 3) ^ _rn 
^ 4 4r2 2 . 

which we write in the form 

[-9,2 + V{r)] fir) = |/(r), f{R) = D/^. (132) 

where 

, (m — 1) (m — 3) , , , 

V{r) = -^^ + 1 ^ DO as r^oo. (133) 

Since 

/ ^dr < 00 and / ^rn^dr < oo for some large Vo- (134) 

Classical results on Schrodinger operators ( see [65] ) allow us to get the asymp- 
totics of /(r) as following: 

/±(r) = Cr-i/2e±X(l + o(l)). (135) 
Now since Uo ^ as r — > oo, we conclude that 

/(r) = /_(r) = Cr-V2e-X(l + o(l)). (136) 

Hence 

Uo{r) =Cr-^e-~{l + o{l)) >0. (137) 

Next, we write 

u{x) = j{x)uo{r) (138) 
Let g e C°°{dB) be strictly positive everywhere and let 

u = Kg. (139) 

Denote by gmm = inf 5 and (/max = sup(/. We obviously have 

5minWo <U< 5maxWo. (140) 

Hence, 
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is bounded. Next, we perform a change in polar coordinates {r,d) by setting 
X = rO. Under this change of coordinates, the operator 



-A+--- + 1 



(142) 



becomes 



m — 1 



a, + - - - + 1 - r-2 



(143) 



where is the Laplace-Beltrami operator on 5"* . Since the operator — AH — - — 



m 



+ 1 is rotationally invariant and dgU takes continuously the value d^g on 

dB, using the fact that each dgU arises as infinitesimal rotation, we conclude 
that for every a, d^u is a solution of the boundary value problem (E) (under 
the change of coordinates) with 



OqU = dgg on dB. 



(144) 



Therefore, 



which implies 



d^j = 0{l), VaeN™. 



(145) 



(146) 



Now we need to control some radial derivative of j. In polar coordinates, we 
have 



m — 1 



+ + l-r-2A, 



Uo{r) = 0. (147) 



Write 



m — 1 



dr + --- + l-r ^Ae 



j{r,0)uo{r)=O. (148) 



Using (129) and the product rule of differentiation, (148) becomes 



drUo f m — 1 



-r ^Aej. 



(149) 



Here 
and 



drUo 



Thus, (149) can be written as 



r 1 
■2+^(;)- 



(150) 
(151) 



-r + 0(-) 
r 



J = 0{r-'). 



(152) 
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Let 

<^(r)=r + 0(i). (153) 

We have 

[dr-f{r)]drj = 0{r-^). (154) 

Let 

F(r) = ^ f{t)dt ~ (155) 

Solving (154) , we get 



(156) 



Since 

F{r) - F{s) ~ - < 2r(r - s) for s > r, (157) 
drj cannot tend to ±oo when r — > c», we conclude that C = and 

POO 

d^j = - e-^W-^^**) [0(s-2)] ds = 0(r-3). (158) 

More generally, since 9^ j is a solution of (157) with right hand side 

- r-^d^ (Aej) = Oir-% (159) 
using the same argument as above with j replaced by dgj, we have 

drd^j = 0(r-3). (160) 

Now differentiating 

[dr-f{r)]drd^j = 0{r-^) (161) 

with respect to r, we get 

[dr-f{r)]d^,d^j = 0{r-% (162) 
using again the same argument as before, we get 

d^^d^j = 0(r-4) (163) 
continuing this way, we finally get 

d^d^j = 0{r-'^-^) k = l,2,... (164) 
Going back to a;— coordinates, we get 

d°'jix) = Oi\x\~^°'^), VaGr",a^O. (165) 
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7 Weighted Estimates for the Decay of Correla- 
tion 



In this section, wc propose to get estimates suitable for obtaining the decay of 
the correlation functions. Wc shall first analyze the case where ^1/ and the source 
term g are compactly supported 



7.1 The compactly supported case. 

We shall assume that $ is given by 

3,2 

^{x) = $a(x) = y + X G R^. (166) 

where 

|a"V«'|<Ca, VaeNl'^l. (167) 

Again g will denote a smooth function on with lattice support Sg = T. We 
shall identify g with g defined on and shall assume that 

\d°'Vg\<Ccc VaeNl^l. (168) 

In addition, we shall momentarily assume that ^ is compactly supported in R"^ 
and g is compactly supported in but these assumptions will be relaxed later 
on. Let M be the diagonal matrix 

where p is a weight function on A satisfying 

^ Pj^l ^ gA if J _ • for sojjje x>0. (169) 
- P(j) - ' ' ^ ' 

Assume also that there exists So S (0, 1) such that 

M-iHess$(a;)M > So (170) 

for every M as above. 
Let 

p(i) = e'«^(^'«s) (171) 

where k is a positive. Define 

/ \ i\2 



Let / be the solution of the equation 



(/W)=0. 



31 



Recall that V/ is a solution of the system 

(- A + V$ • V) V/ + Hess$V/ = Vg in M^. (172) 

Let h = (ti). e 

(V (V$ • V/) ,h} = ifA-^tk + $x Jx.x.ifc) (173) 

= (V/,Hess$ti) + V$- V(V/,ti). (174) 

On the other hand, 

(V(A/),ti) = A(V/,ti). 

We therefore have 

{Vg, ii) = (V$ • V - A) (V/, ii) + (V/, Hess$ti) . (175) 
Because V/(x) — > as |a;| — > oo, we consider a point Xo at which 

VieA / 
is maximal. If M is the diagonal matrix 

M = (Si^pd)) 

we have 

(Vg, Mil) = (V$ • V - A) (V/, Mil) + (V/, Hess$Mti) . (176) 
Now choose 

We need the following lemma. 

Lemma 20 Under the assumptions and notations above, the function 

achieves its maximum value at Xg- 
Proof. Let 

C(x) = (V/(x),Mii) (177) 

and 

i,{x) = \^f{x)\l^. (178) 
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Again by the maximum principle, the function C^{x) achieves its maximum at 
some Xo S M . It is easy to see that critical point for ^(x). Moreover, for 

any a G M^, we have 

(a, Hess7r(xo)a) (179) 
= 2 {a, HessC (a;o)a) + 2^1^ fxixj {xo)fxiXk {xo)p{if j ajUk (180) 

= 2 {a, lIessC{xo)a) + 2 ^ p{if . (a;^) j (181) 

Because (a, Hess7r(xo)a) < 0, we must have {a,'Hess({xo)a) < for any a € 
R^. Thus, local maximum for C(x). Moreover, on one hand, we have 

Cixo) > Cixo) = nixo). (182) 
One the other hand, Cauchy-Schwartz gives 

C{xo) < [7Tixo)f' [7rixo)f^ (183) 
< A^o). (184) 

These last two above inequalities imply 

C{xo) = C(xo) (185) 

and the result follows. ■ 

Now using lemma 21 above, we have 

(V$- V- A)(V/(rCo),Mti) >0. 

This, then implies 

{Vg{xo),Mti) > {Vf{xo),ness^Xo)Mti) 

= {MVf{xo),M-'^lIess^{xo)Mti) 
= (ti,M"^Hess$(a;o)Mti) 

> do\Vf{Xo)\l^. 



Thus 



|V/(x„)||^ < l{MVg{xo),h) 



^\\MVg{xo)\\\Vf{xo)\,^^. 



We have almost proved the following proposition 
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Proposition 21 Let g be a smooth function satisfying 

\d'^Vg\ <Cc Va e N'^I (186) 

and $ is as above. If f is the unique CP°— solution of the equation 

-A/ + V$ • V/ = <? - (5) 
(/W)=0, 

then 

^ fl^ {x)e'^'"^^^'^^^ <C Va; e ffi^ 

C and K are positive constants that could possibly depend on the size of the 
support of g but do not depend on A and f. 

Proof. If 

\'^fM\2,p = (187) 
there is nothing to prove otherwise we have 



< 



1/2 



and the result follows. 



Corollary 22 Let g and h be smooth functions on ]R^,anc! MF where F and 
r' ^ A with r n r' = denote respectively the support of g and h and assume 
that g and h satisfy (4.22) . Then under the assumptions of proposition 2, we 
have 

|cov(5, h) I < Ce-'"^^^'^'^^^ (188) 

where C and k are positive constants that do not depend on A, but possibly 
dependent on the size of the supports of g and h. 
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Proof. Using the formula for the representation of the covariance, we have 



\cov{g,h)\ = (^Ai~ Vg-Vh'^ 
= l(V/-V/i)| 



< 



vieA 



< 



ieA 

1/2 



1/2 



1/2 



1/2 



hlM)dn{x) 



-Kd{Sh,Sg) 



Remark 23 This is the higher dimensional version of theorem 1.4 in [8]. No- 
tice that our proof does not require the assumptions (1.17) and (1-19) namely 



|Hess$(x)||£(, ) < C 



and 



||Hess*(a;)||£(;^) < p < 1 
for all p as above. However, we required that $ satisfies 

M-iHess$(a;)M > So 
for some So € (0, 1) and M as above. 

Notice also that the proof does not require any approximation of mean-field type. 

7.2 Relaxing the Compact Support Assumptions. 

We propose now to relax the assumptions of compact support made previously 
on and g. As before, let M be the diagonal matrix 



where p is given by 
and 



M-iHess$(a;)M > 6o 



(189) 
(190) 
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for every M as above. Next, we propose to generalize the results in propositions 
22 without the assumptions of compact support on and g by means of a family 
of cutoff functions. Let us introduce as in [8] a family cutoff functions 

X-Xe (191) 

(£ G [0, 1]) in C^(M) with value in [0, 1] such that 

X = 1 for |t| < 

\x^''\t)\<Ck^ forfceN. 

We could take for instance 

xAt)=f{e\n\t\) 
for a suitable /. We then introduce 

f = Xe(kl)*, a^eK^ (192) 



and 

Recall that 
which implies 
where 



g,{x) = Xe{\x\)9 X e K^. (193) 

-A/ + V$-V/ = 5-<5>t,A. (194) 

(- A + V$ • V) V + Hess^v = Vg (195) 
v = V/. 



Under the transformations 

V = e~*/^u and q = e~*^^ Vgf 



we have 



- A+i-^-^ — 1 O lu + Hess^u = q in W^. (196) 



We first verify that the assumptions on \1/ and g are satisfied by ^'e(a:) and 
ge{x). Namely 

|a"V*|<Ca, VaeNl^l, (197) 
\d°'Vg\<Cc,, VaeNl^^l, (198) 

and 

M^^Hess^M >(5>0, 0<(5<1 (199) 
M shall still denote the diagonal matrix 

M = ('^i,p(i))i,,eA 
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where p is a weight function on satisfying 



^ pW ^ gA [fi^j for some A > 0. (200) 

Using 

Hess* >S-1, 

we obtain immediately 

M'^Hess*e(3j)Af >iS-l) Xeil^l) - Ce (201) 
for all £ and some constant C. Indeed, we know that 

M-iHess*(x)M > ((5- 1). 
For simplicity we shall write 

Xe = X and r = \x\ 
*e(a;) = x(0*(a;) 



(^^.■-^)x'(r)* + 3^^X»* 



p(i) r p(i) V r2 ; ' ^ p(i) r2 

+ 4^^x'(.)*..+4T^(r)*..., 
p{i) r p{i) 



Let a e 



(M~^Hess*£(a;)Mo, a) 



p(0 

> -2- \x'mx)\ - \x"mx)\ - C \x'{r)\ + {6 - 1) x{r)a^ 



r 

2 



> [{S-l)x{r)-sC]a\ 
We conclude that 

M-^Hess*e(a;)M > (5 - l)x(r) - eC 

for all £ > 0. 
It follows that 

M-iHess$e(a;)M > 5 - C£. (202) 
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Now with 6 replaced hy S' = 5 — Ce, we see that 

M-iHess$e(a;)M > < 5' < 1 (203) 

for e small enough. (Notice that e is possibly A— depend)It remains to check 
the assumptions on and To sec that 

|5"V5e| <C + Oa,A(£), Va G nI^I, (204) 

we have 

9e{x)=Xe{r)9{x), xeR^. 
Again let \a\ > 1, using Leibniz's formula, we have 



/3<a 

= \d'^g\ + \gd-^Xe{r)\ + E f^S) I^^XeW^^-^^l • 

With the assumption g{0) = 0, we write 

\9{x)\ < [ Y.\xjg,.{sx)\ds (205) 

< Cgr (207) 
again using the fact that 

ra«x.(r) = 0„{s), 

we get 

\gd"xAr)\ = OaA^). (208) 

observe also that 

(209) 



/3<a 

it then immediately follows from the assumption on g that 

\d"Vg,\ < C„ + Oa,A(£), Va e NI^L (210) 
Similarly, one can prove that 

|a«V*e| < a + Oa,A(£), Va e NI^I, (211) 

Thus and are compactly supported and satisfy all the conditions that 
were previously required on ^' and g. If denotes the family of solutions 
corresponding to the family of data $£ and ge, one can see that converges 
to u in C°° . The proof which based on regularity estimates is given in detail in 
[8], Consequently, the family of solution = e*=U£ converges to v in C°° 
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Proposition 24 lfg{0) = 0, then Proposition 22 holds without the assumptions 
of compact support on \1/ and g. 



Proof. Using proposition 2 we have 



1/2 



d{i,S,) 



<c\s. 



,1/2 



OA{e) Va;e 



The result follows by taking the limit as e — > ■ 
Corollary 25 If g = Xi and h = Xj we get 

\cor{i,j)\ <Ce-^'^^''^^ 
Which shows that we are away from a critical point. 

8 The d- dimensional Kac Model 

An example of a non-quadratic model satisfying the assumptions above is given 
by 



^a{x) = — 2 In cosh 



The summation is over all nearest neighbor sites. 



^(a;) = -2 ^ In cosh 



{X'l ~\~ Xj^ 



i^Xi ~\~ Xj ) 



with V > small enough. 



It then follows that 



yfsinh {^{x^ + Xj\ 
cosh {xi + Xj)\ 

1 



cosh^ [y^ {xi-\-Xj)\ 



cosh^ [y^ {xi + Xk) 



if k = i 

if /e ~ i 
otherwise. 



and 



l*xj<4dy-, 



|*xiXfc| if fc ~ i. 
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Similarly, using the properties of cosh and sinh and the fact that sinh t < cosh t 
for all t one can see that all derivatives of order greater than or equal to one are 
bounded. Now we propose to check that for u small enough, the Kac Hamilto- 
nian satisfies 

M-^B.ess<P{x)M > So 

for some do G (0, 1) and M as above. 
We need the following lemma. 

Lemma 26 (Schur's Lemma- The R and C bound) For each rectangular 
array 

(Cij)l<i<m 
l<j<n 

and each pair of sequence (a;i)i<^<^ and (yj)i<j<„ we have the bound 



< vie ix, 



vi=l 




where R and C are the row sum and column sum maxima defined by 



R = max I 



and 



C = max Yj I 



i=i 



This bound is known as Schur's Lemma, but, ironically, it may be the second 
most famous result with this name. The Schur's decomposition lemma for n x n 
matrices is also known under this name. Nevertheless, this inequality is surely 
the single most commonly used tool for estimating a quadratic form. Going 
back to the example, we have for any a = {ai)i^j^ G K^, 

(M-iHess$Ma,a) 



aiUi 



Now using the Schur's lemma above, we have 

Pii) 



P{3) 



aia-i 



^ E 



aiUi 



< 



/RCa^ 
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where 



and 



R = max 



C = max 
To estimate i?, observe that 



E 



p(i) 



Now using the fact that 

we have 
Hence 

Similarly, we have 
Thus, 



e"'" < 44 < e'" if i ~ j 



E 



< 2du + 2due'^. 



pU) 

R < 2du {1 + e'^) . 
C < 2diy {1 + e") . 



(M-^Hess$Ma,a) > [{1 - 2diy) - 2dv {1 + e"^)]: 



= l-2dv{2- e"-) . 



The result follows by choosing < k < In 2 and v < 



1 



2d(2- e«)' 
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